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^ • We examine a simple, fuel-air, model of combustion in spark ignition (si) engine with 

x/^ . indirect injection. In our two fluid model, variations of fuel mass burned in cycle sequences 

O [ appear due to stochastic fluctuations of a fuel feed amount. We have shown that a small 

[ amphtude of these fluctuations affects considerably the stability of a combustion process 

' strongly depending on the quality of air- fuel mixture. The largest influence was found 

^ ■ in the limit of a lean combustion. The possible effect of nonlinearities in the combustion 
process has been also discussed. 
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: 1 Introduction 

Cyclic combustion variability, found in 19th century by Clerk (1886) in all spark ignition 
(si) engines, has attracted a great interest of researchers during last years (Heywood 1988, 
Daily 1988, Foakes et al. 1993, Chew et al. 1994, Hu 1996, Daw et al. 1996, 1998, 2000, 
Letellier et al. 1997, Green Jr. et al. 1999, Rocha-Martinez et al. 2002, Wendeker et al. 
2003, 2004, Kamihski et al. 2004). Its elimination would give 10% increase in the power 
output of the engine. The main sources of cyclic variability were classifled by Heywood 
(1988) as the aerodynamics in the cylinder during combustion, the amount of fuel, air 
and recycled exhaust gases supplied to the cylinder and a mixture composition near the 
spark plug. 

The key source of their existence may be associated with either stochastic disturbances 
(Roberts et al. 1997, Wendeker et al. 2000) or nonlinear dynamics (Daw et al 1996, 1998) 
of the combustion process. Daw et al. (1996, 1998) and more recently Wendeker et al. 
(2003, 2004) have done the nonlinear analysis of the experimental data of such a process. 
Various attempts have been done to explain this phenomenon Shen et al. (1996) mod- 



Table 1: Constants and variables of the model. 
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elled a kernel and a front of flame in the cylinder. Their motions and interactions with 
cylinder chamber walls influenced the region of combustion leading to turbulent behaviour. 
Hu (1996) developed a phenomenological combustion model and examined the system an- 
swer on small variations of different combustion parameters as changes in fuel-air mixture. 
Stochastic models of internal combustion basing on residual gases where considered by 
Daw et al. (1996, 1998, 2000), Roberts et al. (1997), Green Jr. et al. (1999), and re- 
cently by Rocha-Martinez et al. (2002), who examined multi- input combustion models. 
Daw and collaborators described the combustion process using a recurrence model with 
nonlinear combustion efficiency and invented symbolic analysis for combustion stability. 

In this paper follow these investigations with a new model assuming that variations of 
a fresh fuel amount is the most common source of instability in indirect injection. 

The present paper is organised as follows. After the introduction in the present section 
we define the model by a set of difference equations in the next section (Sec. 2). This 
model, in deterministic and stochastic forms, will be applied in Sec. 3, where we analyse 
the oscillations of burned mass. Finally we derive conclusions and remarks in Sec. 4. 

2 Two fluid model of fuel-air mixture combustion 

Starting from fuel-air mixture we define the time evolution of the corresponding amounts. 
Namely, we will follow the time histories of the masses of fuel m/, and air rua- 




Figure 1: The combustion curve Am(A) for the constant fresh air feed 6nia = 200 mg. 



Firstly, we assume the initial value of ma{i), iTifii) and automatically their ratio r{i): 

... rngji) 
mf{i) 

for i — 0. 

Secondly, depending on parameter r with reference to a stoichiometric constant s we 
have two possible cases: fuel and air deficit, respectively. For a deterministic model, the 
first case lead to 

r(i) > 1/s (2) 
we calculate next masses using following difference equations: 

mf{i + l) = a ^m/(i) ma{i)^ + Srrif 

ma{i + l) = Sruf, (3) 

where a is the residual gas fraction of the engine, 5m f and Snia denotes fresh fuel and air 
amounts added in each combustion cycle i. In the opposite (to Eq. 2) case 

r{i) < 1/s (4) 

we use the different formula 

m/(?' + 1) = SrUf 

ma{i + l) = a{ma{i) - smf{i)) + 5ma (5) 

Note that variables nia and m/ are the minimal set of our interest. From the above 
equations one can easily calculate other interesting quantities as the combusted masses of 
fuel Am/(i) and air Ama{i) and air-fuel equivalence ratio before each combustion event 

mf{z) + aAmf{i — 1) 
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Figure 2: The dependence of A on sequential cycles i for deterministic (a) and stochastic 
(b) processes. The fresh air amount was assumed to be SrUa = 200 mg while the fresh 
fuel amount varies: 6mf = 13.50 mg , 14.63 mg, 21.00 mg starting from the top curve, 
respectively. 



Basing on experimental results we use the additional necessary condition (Kowalewicz 
1984) of combustion process 

0.6 < A < 1.3. (7) 

For better clarity our notations of system parameters: constants and variables are sum- 
marised in Tab. 1. 

Basing on the relations (Eqs. 1-7) we plotted the combustion curve for the assumed con- 
stant fresh air feed SrUa = 200 mg. This value will be used for all simulations throughout 
this paper. 

Finally, in the case of stochastic injection, instead of constant 5mf{i) (Eqs. 3 and 5) 
(for each cycle i) we introduce its mean value Srufo = const., while 6m f in the following 
way: 

Smf{i) = Smfo + -fN{0,l,i), (8) 

where A^(0, 1, i) represents random number generator giving a sequence i of numbers with 
a unit-standard deviation of normal (Gaussian) distribution and the nodal mean. The 
scaling factor 7 = cxm/ corresponds to the mean standard deviation of the fuel injection 
amount. The cyclic variation of 6mf{i) can be associated with such phenomena as fuel 
vaporisation and fuel-injector variations. 

3 Oscillations of burned fuel mass 

Here we describe the results of simulations. Using Eqs. 1-8 we have performed recursive 
calculations for deterministic and stochastic conditions and obtained time histories of 
various system parameters: mj, rria, Arrif, Am^ and A. The results for A are shown in 
Fig. 2. The upper panel (Fig. 2a), corresponding to deterministic combustion for three 
different values of fuel injection parameter 6mf, shows A as straight lines versus cycle i, 
while the lower (Fig. 2b) one reflects the variations of A in stochastic conditions. The 
order of curves appearing in the Fig. 2b is the same as in Fig. 2a stating from the 
smallest value of considered fuel injection amounts from the top. To get a more clear 
insight of random fuel injection on the engine dynamics , in our stochastic simulations, 
we assumed standard deviation of its mean value 7 = amf = 0.1 Snifo to be high enough. 
In following calculations it was equal to 10%. The obtained results clearly indicate that 
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Figure 3: The dependence of burned fuel mass on sequential cycles i for deterministic (a) 
and stochastic (b-d) processes. SrUa = 200 mg while Snif takes different values: 13.50 mg 
, 14.63 mg, 21.00 mg denoted in particular figures. 



the fluctuations of A are growing with larger A. This can be also found by analytical 
evaluation of Eq. 6. It is not difficult to check that 

cr\ ~ AV^/. (9) 

The results for burned fuel mass Am/ are presented in Fig. 3. Starting from determin- 
istic conditions {Smf = 5m fo = const.) wc obtain the constant fraction of the burned fuel 
mass Amf represented by the three straight lines in Fig. 3a. All lines are lying very close 
to each other and they arc hardly distinguishable in Fig. 3a, but actual numbers shows 
clearly that Am/ increases slightly with growing 5m/ (Am/ = 13.50 mg for 6mf = 13.50 
mg while Am/ = 13.67 mg for 5mf — 14.63 mg and 21.00 mg) account for combustion 
constraints (Eqs. 1-7) and combustion curve (Fig. 1). 

In Fig 3 b-c we show the same. Am/, for the considered case of assumed fuel injection 
{Smfo = 13.50 mg - Fig. 3b, 6mfo = 14.63 mg - Fig. 3c, 5m/o = 21.00 mg - Fig. 3d) 
and stochastic conditions. Due to different magnitudes parameter A fluctuations, and 
dependence of combustion curve Fig. 1 it is not surprise that the fluctuations of Am/ 
have different character in all these cases. For lean combustion, which is a stable process 
in deterministic case, the fuel injection fluctuations introduce considerable instabilities to 
the combustion process leading to the suppression of combustion because in some cycles 
(Fig. 3b) where A is larger that 1.3. 

Then Equation 7 is not satisfied. In the next case (Fig. 3c) the effect of stochasticity 
is much smaller. Here we have optimal air-fuel mixture. First of all one should note 
that fluctuations of A are smaller than in previous case (Fig. 2b). Moreover A oscillate 
around the region (A ~ 1) in combustion curve (Fig. 1) which does not have big changes 
comparing to previous case. Finally, Fig. 3d shows the sequence of Am/ for the large 
5m fo (rich fuel-air mixture) . The fluctuations of A are the smallest of all three ones but 
A ~ 0.7 causes suppressions of combustions in some cycles similarly to the case shown in 
Fig. 3b. 
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Figure 4: Return maps for burned fuel mass: mj(i + 1) versus mi{i), where i denotes a 
present cycle for stochastic process (7 = 10%). SrUa — 200 mg while Sruf takes different 
values: 13.50 mg , 14.63 mg, 21.00 mg denoted in particular figures. Note, arrows in Figs. 
4c and f indicate singular points on return maps for Sruf— 21.00 mg. 




Figure 5: Bifurcation diagrams. Snif against noise (Fig. 5a-c) and against the fresh fuel 
amount constant 5m/ for a — 0.08. Arrows indicate a bifurcation regarding to misfires 
appearance. 
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Figure 6: Bifurcation diagrams: Am/ against a residual gas fraction a for various fresh 
fuel amounts 5m f. 



In Figs. 4a-c we have plotted return maps of all considered stochastic cases Figs. 3b-d. 
Additionally, we have shown similar plots for higher residual gas fraction a = 0.16 (Figs. 
4d-f). Note that differences between both cases {a =0.08 and 0.16) are small. It is also 
worth to note this results are agreement with diagrams of experimental heat release data 
for engines Kohler and Quad4 presented by Green Jr. and coworkers (Fig. 1 in Green ct 
al. 1999) for three values of an equivalence ration. Of course heat release can be treated, 
in the first approximation, as proportional to burned mass through a heating constant of 
fuel. The only important difference visible between our simulated and experimental data 
is that the second is more smeared. This could be connected with our model assumptions 
basing on two components. In reality thereare more components (Heywood 1988, Rocha- 
Martinez et al. 2002) and consequently instead of the simple relation between burned mass 
and heat release there is more sophisticated one. The second, more important, reason may 
lie in our assumed combustion characteristics (Fig. 1). It eliminates partial combustion 
in the initial stage and slow development of flame for smaller A may be important for 
combustion process dynamics (Hu 1997). Our results are also close to that obtained from 
simulations by Daw et al. (Daw et al. 2000, Fig. 9). 

For better clarity we have also plotted bifurcation diagrams with respect to stochastic 
parameter 7 e [0%, 22%] (Figs. 5a-c) as well as a fresh fuel feeding constant Smf E [8, 22] 
(Figs. 5d) and 7 = 10%. Note, single points in vertical cross sections indicate stable 
combustion, multiple points combustion instability with misfire and possible oscillations. 
Finally, dark regions identify combustion with stochastic oscillations and possible inter- 
mittent misfire. 

In Figs. 6a-c we have shown bifurcation diagrams against residual gas fraction a 
{a G [0.05,0.25]). It is clear that this parameter has a small infiuence on combustion 
fluctuations. We have not observed any qualitative changes for all considered fresh fuel 
amounts 6m f. 



4 Conclusions 



In this paper we examined the origin of burned mass fluctuations in a simple model of 
combustion. In case of stochastic conditions we have shown that depending on the quality 
of fuel- air mixture the final effect is different. The worse situation is for lean combustion. 
The consequences of it can be observed for idle speed regime of engine work. Unstable 
engine work, interrupted by the cycles with misfire, leads to a large increase of fuel use. 

Although the presented two component model is very simple it can reflect the underlying 
nature of engine working conditions. In spite of fact that the model is characterised by 
the nonlinear transform (Eqs. 3, 5 and 7) similar to logistic one, we have not found any 
chaotic region. On the other hand we do find qualitative and quantitative features given 
in the earlier experimental works (Green Jr. et al. 1999). However we cannot exclude 
that chaotic solutions can be found for other engine parameters. 

The other strong limitation was concerned with the sharp edges of combustion curve 
(Am/ versus A in Fig.l) modefled by a sharp decay (step function). We used such an 
approximation simplest one 

It can be improved but modelling with the exponential growth exp(— l/x). Such as- 
sumption would be more realistic as it would correspond to partial combustion where the 
mixture of air and fuel is changing its properties inside the cylinder. Prom a physical point 
of view mixture gasoline-air is not uniform before ignition and that can cause nonuniform 
combustion smearing the edges of the combustion curve (Fig. 1). In such a case slow 
development of early flame can also vary from time to time. Going in that direction one 
can also include additional dimensions to incorporate diffusion phenomenon (Abel et al. 
2001). 

However assumptions about the exponential dependence of combustion curve may in- 
volve an averaged effect. We think that it led to chaotic behaviour in earlier papers (Daw 
et al. 1996, 1998, Green Jr. et al. 1999). Calculations considering this effect in our model 
are in progress and the results will be reported in a separate future publication. 
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